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Abstract 



We studied oscillatory behavior of critical amplitudes for the Gaussian 

model on a hierarchical structure presented by a modified Sierpinski gasket 

lattice. This model is known to display non-standard critical behavior on the 

lattice under study. The leading singular behavior of the correlation length ^ 

near the critical coupling K = Kc is modulated by a function which is periodic 

in In I ln(iCc — K)\. We have also shown that the common finite-size scaling 

hypothesis, according to which for a finite system at criticality ^ should be of 

the order of the size of system, is not applicable in this case. As a consequence 

of this, the exact form of the leading singular behavior of ^ differs from the 

one described earlier (which was based on the finite-size scaling assumption). 
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I. INTRODUCTION 

In the near past a considerable research activity has been devoted to the studies of 
recursion relations which have a singular structure near the pertinent fixed points [1-6]. It 
was found that, under certain conditions, these singularities can lead to an unusual critical 
behavior of relevant physical quantities. In particular, it has been shown that the mean 
end-to-end distance R^ of a simple ideal polymer chain on some hierarchical structures can 
grow more slowly than any power of its length A^. This effect has been termed localization, 
and it has been attributed to an entropic trapping of the polymer chain []T|: In order to 
maximize the entropy, it is advantageous for a chain to visit the lattice sites of the highest 
coordination number preferentially. These sites act, therefore, as entropic traps preventing 
the swelling of the chain. 

It is well known that statistics of an ideal polymer chain on lattices can be captured by 
means of a suitable Gaussian model. Using this connection a number of interesting results 
for the polymer model have been derived by studying the singular structure of associated 
recursion relations for the Gaussian model. Let us call here just one example: It has been 
argued, by using finite-size scaling arguments and an analytical study of a pertinent mapping, 
that the mean end-to-end distance of an ideal chain on the modified Sierpinski gasket (see 
Fig. 1) follows the logarithmic asymptotic law: i?^ ~ In A^, with $ = In 3/ In 2. As it has 
been emphasized |^, however, it was very difficult to check this result numerically, due to 
the oscillations in the values of the Gaussian correlation length ^. In fact, these oscillations 
are so pronounced that they mask even leading singular behavior of ^. 

This type of corrections to the leading asymptotic behavior near criticality were studied 
quite early p| , within the framework of the renormalization group approach. Recently, they 
were observed in many systems displaying the power law singularities, and they were related 
to the concept of discrete scale invariance 0. This motivated us to reconsider here the 
Gaussian model on a modified SG lattice. In this case the hierarchical structure of the 
lattice leads to the oscillations of critical amplitudes of various quantities close to criticality. 



If one presents them on an appropriate scale, it turns out that they become regularly spaced. 

We have found that the above mentioned oscillations are universal, and that they can be 
described in terms of some simple functions which are periodic in the variable In | ln{6K) \ 
(with 6K = Kc — K being the distance from the critical point). This is in contrast to the 
critical behavior of usual systems (displaying the power law singularities), in which case 
corresponding variable has the form ln{6K). What is perhaps even more interesting, we 
have found that our numerical values of the correlation length do not fit the form ^o In i^^)^ 
with the above reported value $ = In 3/ In 2. We have shown instead, using an asymptotic 
matching, that one has to take $ = ln(3/2)/ In 2 in order to have a proper agreement between 
analytical and numerical results. Consequently, the common finite-size scaling arguments, 
which were used in the previous studies of this model [||,|],0, are not applicable in this case. 

In Sec. II we present our model, recall some previously obtained results, and examine 
oscillatory behavior of some critical amplitudes. Our conclusions are given in Sec III. 

II. THE MODEL AND ITS ANALYSIS 

We consider here the usual zero-field Gaussian model on a hierarchical structure which 
is presented by a modified Sierpinski gasket (SG) lattice of base 6 = 3 (Fig. 1). Partition 
function of this model has a simple form, 

OO -I 

Z{K) = J2dS,... dSN exp[-- J2S^ + kY: S,S,], (1) 

—CO i <ij> 

where Si is the continuous spin variable at site i, K = J/ksT where T is temperature and 
J represents interaction between each nearest-neighbor pair of Gaussian spins, while < ij > 
denotes the summation over all such pairs. As we have shown in a previous paper ||^, the 
rth order partition function can be expressed in terms of three parameters A^'^\ B^^\ and 
/}('') which obey a set of recursion relations and suitable initial conditions. These relations 
are somewhat cumbersome (see Eqs. (34)- (35) of [^]), and we will not repeat them here. 
Let us note, nevertheless, that they have a singular structure near the relevant fixed point. 



which does not allow us to make a common fixed-point analysis. As detailed in [^], for the 

critical value K^. = 0.227 148 ... of the interaction strength K, all successive iterations of 

A^"^^ and B^^^ lie on an invariant line, which starts at the point {A^^^ = 0, B^^^ = Kc) and 

ends at the fixed point [A* = 1/6,5* = 0). An asymptotic equation of this line, which is 

valid near the fixed point, has been found perturbativelly |^. One can show that along this 

line {6K = 0) parameter B renormalizes according to the law 

21 
B' = 2V3B^ + (18 + 12V3)B^ - -jVSB^ + ..., (2) 

while away from it {6K > 0) this parameter follows the law: B' ~ B^. It is evident that 
the solution Bj. ~ /i'^' (where /i = fi{6K) < 1) well fits in with the latter condition. On the 
other hand, an asymptotic solution of (2) can be expressed as a power series in k^'^, 

r(0 = ^^2'- _ 2 + V3 2-T- , 312 + 193^3 3.2. 

6 8 384 ' ^ ^ 

where < k < 1 is a constant that can be determined numerically. Although this solution 

has been derived for 6K = 0, it holds also for a finite but very small value of 6K (0 < 

6K ^1), provided r ^ Tq, where Tq ^ 1 is the number of iterations one can make along 

the invariant line before going away from it. These two regimes are separated by a rather 

narrow crossover region, which makes it possible to apply the asymptotic matching: 

H^^° = exp[ln(/i) 3"°] ~ k^^° = exp[ln(fi:) 2"°], (4) 

i.e. 

i-{D\ (5) 

where ^{SK) = —1/ ln[/x(5i^)] stands for the correlation length of the model 0. This finding 
is in contrast to the usual finite-size scaling expectation, according to which the correlation 
length of a finite system at criticality should be of the order of system's size [^ ~ 3'"°). 

The number of iteration vq along the invariant line depends on the value of 6K and 
can be estimated from the obvious relation: B^'^''\SK) ^ B^^°\0) + '^^1? 5K. Indeed, 



SK=a 



taking into account (3) and relation dB^'^^^dK ~ 2^° (see 0), we find: 



K,^ ° ~ 2'°6K, or 2'^" ln(fi:) ~ ln{5K). (6) 

This, together with (6), leads to a logarithmic singular behavior of ^ 

e~|ln(<5ir)|^ with <^ = ^^^- (7) 

This differs (in the value of $) from the earlier reported results which have been de- 
rived by using the finite-size scaling assumption: ^ ~ 3^°. This difference is caused by the 
peculiar behavior of the parameter B and correlation function - in our case they decrease 
exponentially to zero even at criticality (see (3) and note 0). 

In order to provide some further insight into the critical behavior of the model, we will 
also examine it numerically. Thus, using a huge precision, we have been able to come very 
close to the critical point {6K/Kc < 10^'^°°°). This allows us to calculate the correlation 
function and associated correlation length in a wide region (on a logarithmic scale) around 
the critical point [|10[- Our results are presented in Fig. 2(a), where the scaled correlation 



length (critical amplitude ^o) ^o = 'CI ln(5i^)|~* is displayed as a function of In | \n{6K)\. The 
overall behavior of ^o is highly sensitive to the precise value of $. For example, average value 
■Co = {(,o,max + ^o,mm)/2 of ,^0 appears to be a constant for sufficiently small values of 6K, 
with the above quoted value of $ (see Fig. 2(a)), while ^q becomes unstable under a small 
change of $. This provides a good criterion for a numerical calculation of $. Indeed, in this 
way we have been able to determine $ with four correct digits, and a further improvement 
depends on the possibility to approach the fixed point still more closely. This should be 
contrasted with the straightforward procedure [Q: A plot of In^ versus In | ln((5-ft')| leads 
to numerical estimates of $ which oscillate with large amplitudes around the exact value, 
independent on the distance 6K from the critical point. 

It seems that ^o represents a simple periodic function of In | ln((5i^)|. Period of this 
function, estimated numerically, is found to be in excellent agreement with the theoretical 
value r = In 2. Perhaps the simplest way to understand this is to adopt the following point 
of view: One can regard (7) as a 'pure' power law, with the scaling variable | ln((5i^)| (rather 
than SK) and a 'critical exponent' $. In the same spirit, one can interpret the relation (6) 



as A*"" ~ I ln((5i^)|, with A = 2 playing the role of a 'thermal' eigenvalue, while relation (5) 
provides an 'effective' spacial scaling ratio 3/2. It is clear then, from standard theory of 
log-periodic corrections to the power law scaling 0, that critical amplitude ,^o should be a 
periodic function in In | ln(5i^)|, with period r = In A = ln2. 

We have also analysed the critical behavior of the first derivative of the free energy 
density with respect to K (internal energy E). Using the approach described in |^, we have 
found that this quantity exhibits an interesting confluent singularity, 

i^^_L|in(5ir)|-*, with ^ = ||^, (8) 

which corresponds to a first order phase transition. As in the case of the correlation length, 
we have studied the internal energy numerically. Our results are displayed in Fig. 2(b), 
where we presented the scaled energy (i.e. critical amplitude Eq) Eq = 6K\\iii{6K)\'^E 
as a function of In | ln((5ii')|. It is evident that Eq represents a simple periodic function 
of ln|ln(5i^)|, while its period is in agreement with the above quoted theoretical value 
(r = In 2). At the same time this analysis provides a good numerical check of the form (8) 
of the energy leading singularity. 

III. CONCLUSION 

In this paper we have studied critical behavior of the Gaussian model on a modified 
SG. We have shown that both correlation length and energy critical amplitudes exhibit very 
pronounced oscillations near the critical coupling Kc- These oscillations can be described 
in terms of some simple functions which are periodic in In | ln((5i^)|. Period r of these 
functions is found to be determined by a universal quantity which governs critical behavior 
of the model. Knowledge of these functions is very useful because it provides more precise 
description of the critical behavior of quantities under consideration. 

Having determined the basic properties of these functions, we have been able to make 
a precise numerical check of the exact form of the leading singular behavior of ^ and E. 



In particular, this allowed us to notice an inaccuracy in a previously described leading 
asymptotic form of C, near criticality |1|,0J^. It seems that this discrepancy stems from the 
inapplicability of the standard finite-size scaling assumption in this model. Indeed, using 
a simple technique, not relying on this assumption (an asymptotic matching), we derive 
somewhat different singular behavior of ^ (see (7)), which turns out to be in excellent 
agreement with the acquired numerical findings (Fig. 2). This example points out that an 
'ad hoc' use of finite-size scaling assumptions could be questionable sometimes, and that 
one has to use them with caution in general case. 
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FIGURES 

FIG. 1. First two stages in the iterative construction of the modified 6 = 3 Sierpinski gasket. 

FIG. 2. (a) The correlation length critical amplitude .^O) in units of lattice constant, as a 
function of ln\ln{6K)\ (^o = C\M^K)\~'^, <I> = ln(3/2)/ln2). (b) The energy critical amplitude 
Eq, in units of the interaction strength J, as a function of the same variable {Eq = 6K\ ln{6K)\^ E, 
^ = In 6/ In 2). These functions are computed from exact representations of related quantities by 
using a huge numerical precision [10]. Note that the amplitudes of these oscillations are so large 
that they, in fact, mask the leading asymptotic behavior of ^ and E. 
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